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摘 要 : 太阳 高 能 活动 爆发 与 活动 区 内 的 电流 结构 有 着 密切 的 联系 ,安培 ( Ampere ) 3E 
j=- (V Xx B): 是 测量 活动 区 内 视 向 电流 密度 的 理论 基础 。 由 于 实测 的 矢量 磁场 中 不 可 避 
免 地 存在 随机 噪声 ， 因此， 应 用 安培 定律 的 不 同形 式 计算 的 电流 密度 存在 显著 的 差异 。 为 
了 比较 不 同形 式 计算 结果 的 差异 并 从 中 探索 一 种 实用 的 电流 计算 方法 ， 基 于 太阳 动力 学 天 
文 台 (Solar Dynamic Observatory, SDO)/HMI(Helioseismic and Magnetic Imager)1£201 1 
年 2 月 15 日 测量 的 活动 区 AR11158 的 矢量 磁 图 ,运用 安培 定律 的 微分 算法 和 积分 算法 分 别 
计算 了 活动 区 内 视 向 电流 密度 的 分 布 图 。 结果 显示 ,微分 算法 获得 的 视 向 电流 密度 j, 分 布 图 
受 随 机 噪声 的 影响 要 远 比 积分 算法 获得 的 结果 大 ， 电流 分 布 图 中 的 电流 结构 没有 积分 算法 
获得 的 结果 清晰 。 另 外 ， 在 扩大 积分 环 路 半径 的 情况 下 ， 所 计算 的 电流 分 布 图 里 的 噪声 信 
号 快速 降低 ， 从 而 使 计算 的 视 向 电流 分 布 图 中 的 电流 结构 更 清晰 。 但 是 当 继 续 扩大 积分 环 
路 半 笃 时， 在 获得 清晰 电流 分 布 图 的 同时 ， 电 流 分 布 图 的 部 分 精细 结构 也 随 之 失真 。 该 研 
究 结果 论证 了 适当 扩大 积分 环 路 计算 视 向 电流 分 布 图 可 以 降低 计算 结果 受 随 机 噪声 的 影响 ， 
从 而 获得 清晰 真实 的 视 向 电流 分 布 图 , 但 是 积分 路 径 的 半径 过 大 在 消除 噪声 影响 的 同时 会 
丢失 电流 分 布 中 的 一 些 精细 结构 。 因 此 在 实际 计算 电流 的 过 程 中 ， 应 该 利用 高 分 辨 率 的 矢 


量 磁 图 ， 选 定 合适 的 积分 路 笃 运 用 安培 定律 的 积分 算法 来 计算 活动 区 的 视 向 电流 ， 从 而 帮 
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助 我们 探索 光斑 爆发 与 活动 区 内 电流 结构 的 关系 。 
关键 词 : 活动 区 ; HM ; 矢量 磁场 ; 电流 密度 ; 安培 定律 


中 图 分 类 号 : P182 ”文献 标识 码 : A ”文章 编号 : 1672-7673 ( 2021 ) 02 


太阳 泡 广 是 太阳 大 气 中 发 生 的 高 能 爆发 现象 ， 其 爆发 过 程 中 常 伴随 日 罗 物 质 抛射 
(Coronal Mass Ejection, CME) 现象 并 向 行 星际 释放 大 量 的 高 能 粒子 流 和 磁 云 ， 成 为 引起 
灾害 性 空间 天 气 的 主要 因素 "”。 因 此 , 研究 炊 斑 爆发 机 制 对 有 效 预防 空间 灾害 性 天 气 有 非常 
重要 的 意义 。 研 究 表 明 ， 磁 流 浮现 和 对 消 ”、 黑 子 旋转 ”*"、 剪 切 运 动 ” 是 触发 大 光斑 爆发 的 
3 种 主要 机 制 ， 在 对 活动 区 AR12673 BABE (包含 一 个 X9. 3 FOOSE 的 研究 中 “提出 了 
一 种 新 的 光斑 爆发 的 触发 机 制 ， 即 由 先前 存在 的 磁场 的 阻挡 引发 的 爆发 模型 。 活 动 区 
AR12673 在 2017 年 9 月 4 日 至 10 日， 产生 了 包括 4 个 X 级 耀 焉 、27 个 M 级 耀 班 和 57 个 C 
ZO PETE PIT] 88 个 耀 班 ， 强 度 最 大 的 两 个 耀 斑 分 别 发 生 在 9 月 6 日 (X9.3 级 ) 和 9 月 10 
H (X8.2 级 )，X9. 3 
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级 次 谢 是 第 24 个 太阳 活动 周 中 最 强 的 耀 班 ， 其 爆发 会 产生 很 强 的 地 球 物理 效应 。 因 此 ， 人 
们 对 该 活动 区 进行 了 大 量 研究 , 相关 研究 结果 表明 , 剪 切 运动 和 黑子 旋转 在 触发 两 个 连续 的 
X 级 兽王 和 日 园 物 质 抛 射 中 具有 重要 作用 ”"， 炮 斑 能 量 释 放 发 生 在 中 性 线 上 方 的 低 剪 切 磁 环 
区 域 ， 光 球 的 增 量 是 由 于 次 评 产 生 的 高 能 粒子 友 击 光 球 产生 的 ”。 在 泡 六 爆发 过 程 中 ， 磁 重 
联 被 认为 是 能 量 快 速 释放 的 有 效 机 制 "， 磁 重 联通 常 发 生 在 剪 切 程度 高 的 磁场 区 域 TER 
重 联 过 程 中 ， 磁 拓扑 结构 发 生 改变 ， 磁 能 转化 为 加 热 等 离子 体 的 热能 和 加 速 粒子 的 动能 。 通 
常 认为 , 位 重 联 最 初 的 能 量 释 放 主 要 发 生 在 日 郡 中, 瀛 斑 的 爆发 主要 由 人 厂 重 联 释放 的 自由 磁 
能 驱动 。 自 由 人 磁 能 是 指 总 磁 能 与 势 场 磁 能 之 间 的 差 值 , 爆发 商议 的 活动 区 磁 位 形 偏离 势 场 越 
大 , 活动 区 的 总 磁 能 与 势 场 磁 能 的 差 值 束 越 大 , 活动 区 的 非 势 性 就 越 强 , 储存 的 自由 磁 能 就 
RE, 炮 广 产生 的 几率 就 越 大 ”。 光 球 及 光 球 下 方 磁 力 线 的 剪 切 和 扭 绞 , 新 浮现 的 磁 流 等 都 
会 形成 活动 区 磁场 的 非 势 性 ”。 由 于 磁 中 性 线 处 的 磁场 梯度 较 大 , 磁场 的 纵向 分 量 无 明显 变 
化 ， 但 水 平 磁场 高 度 剪 切 ， 在 磁 压 力 的 作用 下 产生 大 量 的 自由 磁 能 ， 所 以 多 数 粮 斑 趋向 于 产 
生 在 磁 中 性 线 附 近 ”“。 对 5 个 X 级 耀 斑 的 研究 ”发 现 沿 磁 中 性 线 的 磁场 剪 切 在 耀 斑 之 后 有 
明显 增强 。 最 近 一 些 基 于 高 分 辨 率 观测 获得 的 矢量 人 磁 图 的 研究 表明 , 在 沿 其 主要 磁 中 性 线 有 
强 剪 切 的 活动 区 中 存在 大 量 净 电 流 “。 文 [17] 对 活动 区 电流 分 布 的 研究 表明 , REESE HE SG 
物质 抛射 和 兆 广 的 活动 区 与 生 俱 来 共有 大 量 净 电 流 ， 与 最 近 的 观测 结果 一 致 , Th EEE TH A) 
位 置 和 形态 与 光 球 层 强 电流 带 有 密切 关系 , 在 耀 斑 演化 过 程 中 , 耀 斑 带 的 形态 改变 可 能 与 电 
流 带 相关 ”。 太 阳 活 动 区 中 电流 密度 较 高 的 区 域 ” 与 Hs 兆 斑 核 块 及 具有 极 大 剪 切 的 磁场 区 
域 相 重合 “， 这 个 结果 支持 了 炊 广 是 由 强 电 流 引 起 的 等 离子 体 的 不 稳定 性 触发 的 理论 ”。 
这 些 研究 都 证 实 了 电流 可 能 在 炮 广 爆发 之 前 就 已 存在 , 并 在 光斑 爆发 中 起 重要 的 作用 , 活动 
区 中 电流 带 与 次 闹 带 之 间 密 切 相 关 。 因此, 计算 活动 区 磁 中 性 线 两 侧 强 磁场 中 的 电流 分 布 对 
于 预测 光斑 爆发 的 位 置 和 形态 具有 重要 意义 。 
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随 着 观测 仪器 的 时 间 和 空间 分 辩 率 不 断 提 高 , 人 类 对 太阳 爆发 活动 机 制 有 了 比较 深入 的 


If EP, VERBA E REIS JU ONE 


的 测量 是 基 了 


塞 曼 效 应 进行 的 , 由 于 来 自 


KS 


爆发 活动 的 触发 机 制 非常 重要 。 太 阳 矢 量 磁 场 


日 蛙 、. 过 渡 区 等 太阳 大 和 气 的 偏振 信号 非常 弱 , 因此 ， 


这 些 区 域 的 矢量 磁场 目前 还 很 难 获得 , 目前 唯一 能 够 精确 测量 磁场 的 区 域 只 有 光 球 层 。 基 于 


确 计 算 日 时 中 的 电流 分 布 
层 的 视 向 电流 密度 分 布 。 


o 


测量 的 光 球 层 磁场 结构 ， 人 们 以 此 为 边界 条 件 运 ) 
得 三 维 空间 矢量 人 磁场 的 可 能 结构 。 由 于 无 法 准确 测量 日 蝎 层 的 矢量 侯 场 , 也 就 几乎 不 可 能 准 
基于 精确 测量 的 光 球 矢量 磁 图 , 我 们 可 以 
日 于 观测 仪器 在 测量 矢量 磁场 的 过 程 中 不 可 避免 地 受到 系统 误差 和 


] 非 线性 无 力 场 (NL 


FFF) ”模型 进行 外 推 获 


运用 安培 定律 计算 光 球 


随机 误差 的 影响 ， 其中， 系统 误差 主要 来 源 于 对 太阳 大 气 的 相关 参数 所 做 的 假设 和 简化 、 仪 


器 散光 和 镜 方 的 非 正 常 


引起 ”。 系统 误差 可 以 在 某 种 程度 上 从 原始 数据 


本 


M. RESTI 


F。 随 机 误差 主要 由 大 气 视 宁 度 、 观 测 仪器 电 噪声 等 
消除 或 采用 如 瞳 场 处 理 、 平 场 处 理 等 方法 


减 小 它们 对 结果 的 影响 , 无 法 消除 的 那 部 分 系统 误差 在 整个 视 场 是 连续 可 导 函 数 , 它们 对 电 
流 的 影响 也 有 一 定 规律 。 而 随机 误差 在 不 同 的 空间 格 点 上 随机 变化 ， 


消除 , 通常 服从 正 态 分 布 。 这 给 


它们 不 能 从 原始 数据 中 


我 们 运用 安培 定律 微分 算法 计算 视 向 电流 密度 带 来 很 大 的 障 


但 。 文 [24] 通 过 在 矢量 磁场 中 引入 随机 噪声 , 研究 了 随机 噪声 产生 的 影响 和 空间 分 辨 率 之 间 
的 关系 ， 并 分 析 微 分 法 、 小 环 路 积分 法 和 大 环 路 积分 法 得 到 的 结果 ， 


时 , 沿 一 个 大 环 路 的 安培 定律 的 积分 形式 计算 电流 密度 分 布 得 到 的 结果 最 好 。 为 了 进 
证 梁 红 飞 等 人 的 结果 ， 并 检验 次 斑 带 和 电流 带 的 对 应 关系 ， 我 们 利 月 
月 15 日 观测 到 的 活动 区 AR11158 的 矢量 磁场 ， 运 | 
大 环 路 积分 ) 计算 该 活动 区 的 视 向 电流 分 布 ， 


路 积分 、 中 环 路 积分 、 


也 们 发 现 当 忽略 高 阶 项 
步 验 


H SDO/HMI 在 2011 年 2 


j 安 培 定律 的 4 种 形式 〈 微 分 形式 、 小 环 


将 积分 法 的 积分 环 路 扩 


大 到 适当 的 大 小 进行 比较 。 我 们 还 考察 了 计算 结果 的 电流 密度 极 值 与 光斑 带 的 位 置 和 形态 的 


对 应 关系 ， 以 检验 计算 电流 密度 方法 的 有 效 性 和 合理 性 。 


1 数据 来 源 和 计算 方法 


1.1 数据 来 源 


$ 


SDO (Solar Dynamics Observatory) "'/HMI (Helioseismic and Magnetic Imager) 
2I KARAHAN (Atmospheric Imaging Assembly, AIA) "^ bL A IA iE RAAR CEUV 
Variability Experiment, EVE) 是 搭载 在 SD0 卫星 上 的 3 台 仪 器 。 
太阳 光 球 层 的 矢量 人 磁场 、 光 球 强度 以 及 多 普 勒 速度 等 相关 数据 。2011 年 2 月 15 日 位 于 太阳 


西南 方向 (S20W12) 的 活动 区 AR11158 具 


通过 SDO/HMI 可 以 获得 


有 BY 结构 的 复杂 磁 位 型 , 在 1:44UT-02:06UT 期 间 


该 活动 区 爆发 了 一 个 持续 了 约 22 min 的 X2. 2 级 次席 ,有关 该 活动 区 的 研究 可 参考 文 [28-29]， 


在 本 文中 ， 我 们 通过 SD0/HMI 获得 其 在 炊 广 爆发 期 间 观测 的 磁场 数据 ， 空 间 分 辨 率 为 0.5 


角 秒 /像素 ， 时 间 分 辨 率 为 45 秒 / 帧 。 矢 量 磁场 的 测量 是 基于 塞 曼 效应 运用 侦 振 种 射 转移 方 


程 获得 ， 获 得 矢量 磁场 数据 的 具体 过 程 见 文 
快速 反 演 技术 VFISV(Very Fast 


Inversion 


[30], SDO/HMI 利用 斯 托 克 斯 (Stokes) 矢量 的 
of the Stokes Vector) 5 得 出 矢量 磁 图 。 矢 量 


做 图 的 时 间 分 辨 率 为 12 分 种 / 帧 。 由 于 以 偏振 辐射 转移 方程 为 理论 基础 求解 矢量 磁场 的 过 
太阳 人 磁场 的 测量 带 来 很 大 的 障碍 ， 为 了 解决 这 


程 中 往往 存在 方位 角 


180° 不 确定 性 问 


个 问题 ， 势 场 校准 法 ”、 无 力 场 近似 法 ”” 、 


先后 提出 并 和 


利用 SDO/HMI 
量 磁场 数据 进行 预 处 型 


题 ， 给 


Sb 1 


最 小 能 量 法 ”、 磅 蓓 


导 以 完善 ， 从 而 很 好 地 解决 了 太阳 磁场 测量 中 的 方位 角 180° 不 确定 性 问题 ， 在 


相关 研究 得 以 顺利 进行 。 


[ 测量 的 活动 区 偏振 信号 反 演 矢 引 


EE RH. SDO HI 


s+ [36] > son [37] Ave 
法 ”、 连 续 性 方法 ”等 


队 利 用 能 量 最 小 法 对 矢 


,解决 了 矢量 磁 图 中 的 方位 角 180° 不 确定 性 问题 ”, 使 得 太阳 磁场 的 


chinaXiv:20201 1.00108v1 


ChinaXiv 合 作 期 刊 


SDO/ATA 是 用 于 观测 不 同 高 度 的 太阳 大 和 气 层 的 仪器 ， 可 以 同时 提供 10 个 不 同 波段 的 太 
阳 大 气 信息 。 这 10 个 波段 包括 7 个 极 紫外 波段 CEUV)、2 个 紫外 波段 UV) 和 工 个 可 见 光 


波段 (连续 谱 450 nm). FEHR EUV 波段 包括 
nm, 13. 1 nm; UV 波段 包括 170 nm 和 160 nm. 


30.4nm, 17. 1 nm, 19. 3 nm, 21. 1 nm, 33. 5 nm, 9. 4 
各 波段 在 大 气 中 对 应 的 谱 线 及 观测 区 域 在 文 [38] 


中 进行 了 详细 介绍 。AIA 图 像 的 空间 分 辩 率 为 0.6 角 秒 /像素 , 时 间 分 辨 率 分 别 为 12 秒 / 帧 
CEUV 波段 )，24 $b/üi (UV 波段 )。 本 文 所 用 SDO/AIA 数据 是 170 nm 波段 和 30. 4 nm JE Et, 
其 中 ，170 nm 波段 的 观测 区 域 为 温度 极 小 区 、 光 球 层 ，30. 4 nm 波段 的 观测 区 域 为 色 球 层 、 


过 渡 区 。 


1.2 计算 方法 


1 (a) 为 活动 区 AR11158 Æ 2 H 15 H 02:00:00UT 时 刻 的 纵向 磁 图 ， 其 中 白色 为 正极 


磁场 ， 黑 色 为 负极 磁场 ， 黑 色 框 所 围 


形 负 人 磁极， 下 方 为 一 个 半月 形 的 了 


EWR, ZEW 


也 称 作 纵 向 磁场 的 极 性 变换 线 ， 是 正 负 极 磁场 
02:00:00UT 时 刻 SDO/HMI 观测 到 的 矢量 磁 图 ， 为 了 能 够 清晰 地 看 出 矢量 磁场 的 结构 ， 图 中 


区 域内 存在 一 对 极 性 相反 的 磁极 ,红色 第 头 所 指 位 置 为 
正 负 磁极 分 界 的 磁 中 性 线 如 图 (a) 和 Cb) 中 黄色 虚线 ) 位 置 ， 磁 中 性 线 的 上 方 为 一 个 椭圆 


性 线 附近 磁场 的 梯度 非常 大 。 通 常 磁 中 性 线 


的 分 界线 。 图 1 CO 为 活动 区 在 2 月 15 H 


以 纵向 磁场 为 背景 ,县 加 在 纵向 磁场 上 方 的 小 箭头 表示 横向 磁场 ,其 中 蓝 色 箭头 覆盖 区 域 为 


正极 磁场 ， 红 色 第 头 覆 新 区 域 为 负极 磁场， 第 头 长 度 表示 横向 磁场 的 大 小 ， 第 头 方向 表示 横 
向 磁场 的 方向 ， 灰 度 图 与 图 1 Ca) 完全 相同 ， 均 表示 纵向 磁场 的 大 小 。 从 表示 横向 磁场 的 
箭头 艇 不 难看 出 , 横向 磁 图 具有 旋涡 状 的 结构 , 这 种 结构 意味 着 活动 区 内 可 能 存在 比较 强 的 


Hy. 我 们 利用 这 些 矢 量 磁 图 运用 安培 定律 的 不 同形 式 计 算 活动 区 内 的 电流 密度 分 布 , 并 比 


较 和 分 析 不 同方 法 计算 结果 的 差异 ， 从 而 挑选 一 种 最 优 的 电流 密度 计算 方法 。 太阳 动力 学 天 
文 台 除了 提供 该 活动 区 的 矢量 磁 图 外 ， 还 在 不 同 波段 采集 到 了 活动 区 的 单 色 像 。 1 Cc) 
和 Cd) Jy SDO/AIA 分 别 利用 谱 线 170 nm (01:58: 08UT) 和 30. 4nm (01:58:12UT) 观测 到 
的 2011 年 2 月 15 日 在 活动 区 ARI1158 产生 的 X2. 2 级 耀 斑 的 单 色 像 图 。 从 网 中 可 以 看 出 ， 


不 难 发 现 ， 泡 广 主 要 位 于 人 磁 中 性 线 附 近 。 


泡 议 带 主 要 在 黑色 方 框 区 域内 ， 部 分 结构 位 于 方 框 的 右上 沿 ， 对 比 图 1 Cad, Ceo A Cd) 
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图 1 活动 区 AR11158 的 纵向 磁 图 ， 矢 量 磁 图 和 170 nm PEEL, 30. Anm 波段 的 单 色 像 图 


Fig.1 The Longitudinal magnetogram , vector magnetogram and monochromatic image of the active region 
AR11158 acquired at 170 nm and 30.4 nm wavelength. 
研究 表明 , 在 太阳 活动 区 内 可 能 存在 空间 尺度 相对 较 小 但 物质 密度 很 高 的 片 状 电流 结构 ， 
目前 已 获得 这 些 片 状 电流 结构 大 量 的 观测 数据 , 特别 是 光斑 爆发 期 间 的 观测 数据 , 但 是 电流 
结构 内 的 电流 强度 以 及 具体 的 电流 密度 分 布 特征 目前 还 很 难 准确 测量 到 ””。 我 们 利用 来 自 
HMI 的 光 球 矢量 数据 计算 光 球 表面 的 视 向 电流 密度 产 。 针 对 矢量 磁 图 数据 都 是 离散 的 特点 ， 


通常 情况 下 运用 安培 定律 (Vx B)。 的 差分 形式 


: 1 ,0B OB, 
(uc Ga - 9 (1) 


Ho 


计算 视 向 电流 密度 产 ， 其 中 ，Ho 为 真空 中 的 磁 导 率 。 如 图 2， 利 用 Bx, By 在 最 里 面 的 小 正方 形 


4 条 边 上 的 中 点 AAL y) « BG ytAy). COcAx, y) « DG; 六 ADn 的 值 进行 计算 ， 差分 公式 
变 为 


1 [Cra Ea; Codijn-Bdi, ja 2 
Gali? Lu ee ^ x 


其 中 ，jzo 为 运用 差分 法 计算 的 视 向 电流 密度 。 由 (2) 式 可 以 看 出 ， 空 间 分 辨 率 越 高 ， 空 间 


格 点 宽度 《Ax 和 Ay) 越 小 计算 的 电流 密度 结构 越 精细 。 由 于 实际 所 测 磁场 B: 往往 受 随机 


噪声 影响 而 引起 误差 矢量 $5B， 用 Bt 表示 活动 区 内 磁场 的 真实 值 与 系统 误差 之 和 ， 实 际 所 测 
做 场 可 以 表示 为 Be =Bi+6B, 所 计算 的 视 向 电流 密度 jzo 为 真实 的 视 向 电流 密度 jzo 与 随机 品 
声 产 生 的 误差 电流 密度 6jzo 之 和 ， 即 


, 1_1 (ABy T) 1 (me aes) m 
Jzo m ( Ax Ay 十 Ho N Ax Ay — Jzo F Ojz0 $ (3) 
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HH (3) 式 可 知 ， 当 随机 噪声 5Bx 和 6By 不 变 时 ，Ax 和 Ar 随 着 空间 分 辩 率 的 提高 而 减 小 ， 最 
终 导致 视 向 电流 密度 增加 并 将 真实 电流 密度 淹没 , 所 以 我 们 通常 用 安培 定律 的 积分 形式 
je $iBr /ds (4) 


计算 视 向 电流 密度 ， 其 中 ,Bi 为 实测 横向 磁场 ，Bi = eiBx + ejBy, eiej Bx By WANE x y A 


向 的 单位 矢量 及 实测 横向 磁场 的 两 个 分 量 ; ds 为 积分 路 径 ! 所 围 的 面积 ， 通 常 选择 边 长 为 
7Ax 的 正方 形 为 积分 路 径 ，? 为 空间 格 点 数 , 则 As=z(Axz)。 由 〈4) 式 可 知 ，Ax MERER 
空间 分 辩 率 的 变化 而 变化 ， 分 辨 率 越 高 环 路 需要 扩 得 越 大 ， MM Mn FRE 
通过 分 辨 率 的 高 低 来 选择 适当 大 小 的 积分 环 路 可 以 快速 降低 随机 噪声 对 计算 结果 的 影响 。 通 
常 根 据 分 辨 率 的 大 小 来 定义 小 环 路 、 中 环 路 、 大 环 路 。 


i-3 i-2 i-l i i+l i +2 i43 
一 十 
(x +3Ax J y*3Ay) p 
J y+2Ay) (x+2Ax ,|y+2Ay) p 
jt+l 
C3 
j 
j1 
j-2 
]-3 
y-3Ay) 
图 2 安培 定律 积分 形式 的 积分 路 径 : Ae) TE I APRS EC,» EA Tb SS 2 个 正方 形 为 积分 路 径 Cz， 第 3 个 正方 形 为 


积分 路 径 Cs。 位 于 C1 上 的 4 个 有 星 号 的 点 对 应 于 安培 定律 的 差分 形式 ， 即 基于 横向 磁场 B4 在 ALx+Ax, y]、BLx, y+Ay]、C[x-A 
Xx ,y]、D[x，y-Ayj] 点 的 值 可 通过 Ampere 定律 的 差分 形式 计算 [x，y] 处 的 视 


电流 密度 值 。 


H 


Fig. 2 The integration paths of the Ampere’ s law integral form : The smallest square is the integration path C, from 
the inner to the outside, the second square is the integration path C}, and the third square is the integration 


path C4. The four ‘star’ points which locating in the C, are corresponding to the differential form of Ampere’ s 


law. That is, based on the value of transverse magnetic field B, at A[x+Ax, yl, Bbo y*Ay], C[x-Ax. yl, Dlx, y- 


Ay], the vertical electric current density value at [x, y] can calculate by using the differential form of Ampere' s 


law. 


如 图 2, 环 路 积分 法 是 利用 横向 磁场 Bi 在 积分 路 径 上 的 值 计 算 积分 环 路 所 围 面 积 内 的 平 
均 电 流 密度 ， 得 出 积分 环 路 中 心 0 Cx, 万 点 的 电流 密度 ， 所 得 电流 密度 为 积分 区 域 的 平均 
电流 密度 。 积 分 路 径 沿 逆 时 针 方向 ， 积 分 环 路 以 O Ox, 万 为 中 心 由 小 到 大 (或 由 内 向 外 ) 
依次 称 为 小 环 路 (以 2A x 为 边 长 的 正方 形 ， 积 分 路 径 为 C1 )、 中 环 路 (以 4Az 为 边 长 的 正 
方形 ， 积 分 路 径 为 C2) 和 大 环 路 〈 以 6Ax 为 边 长 的 正方 形 ， 积 分 路 径 为 C3〉 。 其 中 ， 小 环 
路 的 积分 路 径 是 沿 点 (x+Ax ，y+Ay) ， (x-Ax ，y+Ay) , (x-Ax ，y-Ay) 和 (Cx+Ax ， 
y-Ay) 所 在 正方 形 边缘 ， 积 分 路 径 为 C， 对 应 的 视 向 电流 密度 为 zj 。 由 《〈3) 式 可 得 小 环 路 
的 计算 公式 为 


. ub 
(zu 9, Be oa ds 


1 ( i+1 (Bx)n-1,j-1 ix (Eptran-c Bram py 


dnoAxAy “=i 2 = 2 
i+1 (Bx)n-1,j+1 t nj j*1 (By)i-1n-1 *(Gy)i-in 
— yit Enc +P ons ay yit Oii t Olea ay) (5) 


其 中 ， ds-2Ax*2Ay. 


中 环 路 的 积分 路 径 是 沿 点 (x+2Ax,y+2Ay) ， (x-2Ax,y+2Ay) , (x-2Ax, y-2Ay) 和 
(x+2Ax,y-2Ay) 所 在 正方 形 的 边缘 ， 积 分 路 径 为 C2, 对 应 的 视 向 电流 密度 为 zz。 由 〈3) XX 
可 得 中 环 路 的 计算 公式 为 


1 
Gi gu. dI/ds 


1 (xit (Bx)n-1,j-2 +(By)n j-2 Ax + ais (By)i+2n-1+(By)itan Ay 


^ 16uoAxAy “n=i-1 2 n=j-1 2 
i+2 (Bx)n-1,j+2 +(Bx)n,j+2 j+2 (By)i-2n-1 +(By)i-2n 
“Lin=i-1 2 AX. a "— Ay) (6) 


HB, ds=4Ax*4Ay. 


大 环 路 的 积分 路 径 是 沿 点 (x+3Ax,y+3Ay) , (x-3Ax,y+3Ay) , (x-3Ax, y-3Ay) 和 
(x+3Ax, y-3Ay) 所 在 正方 形 的 边缘 ， 积 分 路 径 为 G3, 对 应 的 视 向 电流 密度 为 jz3。 由 (3) 式 
可 得 大 环 路 的 计算 公式 为 


1 
Üza)ij Es $ ,B,.dl/ds 


1 ( i+3 Wie gite Ax + Ye (By), ani +(By) isan Ay 


~ 36uoAxAy “n=i-2 n-j-2 2 
i+3 (Bx)n-1,j+3 +(By)n j+3 j+3 (By)i-3,n-1 +(By)i-3,n 
"Am-i-2 2 Ax — Yn=j-2 > Ay) (7) 


HB, ds-6Ax*6Ay. 


图 3 zs H CREE TRI NE UTA ER BR VV HA FSA Pas cH Tia] C02 : 00 : 00 UT) 
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的 视 向 电流 密度 分 布 图 。 图 中 用 “Dif ”表示 差分 法 ， 用 “Int” 表 示 积 分 法 ,“1C” 表 示 小 
环 路 ,“2C” 表 示 中 环 路 ,“3C” 表 示 大 环 路 。 图 中 白色 和 黑色 分 别 表 示 正 负电 流 密度 分 布 ， 
对 应 的 视 向 电流 密度 值 分 别 为 10. 05A。m-2 和 -0.05 A。m-z。 为 了 便于 比较 和 分 析 ， 我 们 
在 微分 法 和 积分 法 中 所 取 的 视 向 电流 密度 阔 值 相同 ， 均 为 士 0.05 A * m^?. Al 3 Ca)，(b) 
分 别 是 运用 安培 定律 的 微分 法 和 积分 法 《小 环 路 ) 计算 的 视 向 电流 密度 分 布 图 ， 图 3 CoD, 
(d) 是 扩大 积分 环 路 时 计算 的 视 向 电流 密度 分 布 图 。 其 中 〈c) 为 中 环 路 的 计算 结果 ，(d) 
为 大 环 路 的 计算 结果 ,由 于 所 测 矢量 磁 图 中 的 横向 磁场 不 可 避免 地 受 随机 噪声 的 影响 ,微分 
法 计算 的 电流 密度 分 布 ( 如 图 3 (a)) 受 噪声 信号 影响 较 大 , 相 比 之 下 , 小 环 路 的 积分 法 (如 
3 Cb) 计算 的 电流 密度 分 布 受 随 机 噪声 的 影响 相对 较 小 。 当 把 积分 环 路 扩大 到 两 个 (中 
环 路 ， 如 图 3(c)) 时 ， 得 到 的 电流 密度 分 布 图 比 小 环 路 积分 法 计算 的 噪声 信号 更 小 。 当 把 
积分 环 路 扩大 到 三 个 (大 环 路 ， 如 图 3 (d)) 时 ， 虽 然 计 算 的 噪声 信号 更 小 ， 但 电流 密度 分 
布 中 的 部 分 精细 结构 明显 失真 。 因此, 在 考虑 随机 噪声 的 情况 下 ,我 们 选择 扩大 积分 环 路 的 
方法 ， 并 将 积分 环 路 扩大 到 两 个 ， 即 利用 中 环 路 的 积分 方法 计算 视 向 电流 密度 更 好 。 

Mil Dif 2011-02-15 02:00:00 UT © 


int-2011—02—15 02:00:00 UT. 
1C d SUM, 
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图 3 安培 定律 的 微分 算法 和 积分 算法 所 计算 出 的 视 向 电流 密度 分 布 图 . 


Fig.3 The distribution of vertical current density calculated by using the differential algorithm 


and integral algorithm of Ampere”s law. 


由 图 3 Ceo 和 图 1 Ca) 可 以 看 出 ， 强 电流 主要 产生 于 磁 中 性 线 附 近 的 强 磁 场 区 域 ， 与 
强 磁场 有 很 好 的 对 应 关系 。 由 图 1 a), C) 和 Cd) 可 以 看 出 ， 瘤 斑 带 可 能 出 现在 磁 中 性 
线 附 近 的 区 域 ， 由 此 说 明 ， 电 流 带 与 次 斑 带 在 位 置 上 可 能 具有 某 种 对 应 关系 。 如 图 4， 我 们 
通过 做 投影 改正 将 中 环 路 积分 法 (如 图 3(c)) 计算 的 视 向 电流 密度 分 布 与 相同 淋病 时 刻 的 
170nm 波段 和 30. 4nm 波段 的 单 色 像 图 (图 1(c) 和 (d) ) 对 齐 , 图 中 绿色 和 蓝 色 等 高 线 分 别 表 
示 正 负电 流 带 , 对 应 的 视 向 电流 密度 分 别 为 +0. 02 A。m-2 和 -0.02 A * m7. ARER, XE 
斑 带 大 致 出 现在 电流 较 强 的 位 置 ， 此 处 的 自由 磁 能 多 ,释放 的 能 量 也 多 ， 而 且 瘤 斑 带 与 电流 
带 具 有 相似 的 形态 。 
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图 4 电流 带 与 兆 斑 带 的 对 应 关系 , 左 图 为 AIA 170nm 波段 在 2011 年 2 月 15 日 01:58:08UT 时 刻 活动 区 的 单 色 像 图 , A 
图 为 ATA 30. 4nm 波段 在 2011 年 2 H 15 日 01:58:12UT 时 刻 活 动 区 的 单 色 像 图 ， 图 中 的 绿色 等 高 线 和 蓝 色 等 高 线 
分 别 表示 正极 和 负极 电流 带 ， 对 应 的 视 向 电流 密度 值 分 别 为 +0. 02 A。m-2 和 -0.02 A。m-2。 


Fig. 4 The corresponding relationship between the current ribbons and flare ribbons. The left panel is the 


1700A image of the active region observed at 01:58:08UT on February 15, 2011 and the right panel is the 304A 
image of the active region observed at 01:58:12UT on February 15, 2011. The green and blue contour lines 
in the panels represent the positive and negative current values respectively, and the corresponding vertical 


2 


electric current density values are +0.02 A* m ?and -0.02 A* m^? respectively. 


3 讨论 与 结论 


太阳 疙 斑 究 责 在 什么 位 置 爆发 ?其 形态 是 什么 样 的 ? 如 何 寻找 一 种 有 效 预 测 太阳 瀛 斑 
爆发 位 置 和 形态 的 方法 一 直 是 太阳 物理 中 的 一 个 前 沿 课题 。 耀 斑 爆 发 与 活动 区 电流 有 着 密切 
的 联系 , 活动 区 的 电流 主要 表现 为 由 矢量 磁 图 中 横向 磁场 的 变化 产生 的 视 向 电流 。 因 此 ， 选 
择 一 种 有 效 计 算 活动 区 视 向 电流 的 方法 显得 尤为 重要 在 不 考虑 随机 噪声 的 情况 下 , 随 着 空 
间 分 辩 率 的 提高 , 利用 安培 定律 的 微分 形式 和 积分 形式 计算 的 视 向 电流 密度 都 比较 接近 。 在 
考虑 随机 噪声 的 情况 下 ,利用 安培 定律 的 微分 形式 计算 视 向 电流 密度 , 虽然 计算 量 小 , 但 随 
者 空间 分 辨 率 的 提高 《空间 格 点 宽度 Ax 和 Ay 变 小 由 随机 误差 引起 的 电流 密度 的 离散 度 变 
大 , 计算 的 视 向 电流 密度 的 不 确定 性 也 在 增加 一 .即使 将 视 向 电流 密度 分 布 图 进行 平滑 处 理 ， 
在 消除 噪声 信号 的 同时 也 会 把 原来 噪声 信号 覆盖 下 的 电流 信号 消除 。 而 运用 安培 定律 的 积分 
形式 计算 视 向 电流 密度 ,虽然 比 微分 形式 计算 量 大 , 但 可 以 降低 测量 值 的 离散 度 ,使 计算 结 
果 受 随机 误差 的 影响 变 小 ， 若 扩大 积分 环 路 ， 矢 量 磁 场 的 空间 分 辨 率 降 低 〈 空 间 格 点 宽度 
Ax 和 Ay 变 大 )， 可 以 快速 降低 测量 值 的 离散 度 ， 使 万 得 结果 更 准确 。 运 用 安培 定律 积分 形式 
计算 视 向 电流 密度 的 方差 比 微分 形式 计算 的 电流 密度 的 方差 更 小 ”。 当 积分 路 径 超过 所 研究 
矢量 磁 图 的 最 小 空间 分 辩 率 结构 的 线性 大 小 时 , 所 得 电流 密度 分 布 的 部 分 精细 结构 就 会 失真 。 
基于 高 分 辩 率 矢量 人 磁 图 寻找 一 种 合适 的 积分 路 径 计算 视 向 电流 密度 , 能 在 去 除 噪声 信号 的 同 
时 保留 原 有 的 电流 信号 ， 所 得 结果 更 准确 。 

基于 此 理论 ,本 文通 过 比较 计算 视 向 电流 密度 的 几 种 方法 ， 即 安培 定律 的 微分 算法 、 小 
环 路 积分 法 、 中 环 路 积分 法 和 大 环 路 积分 法 〈 如 图 3)， 最 终 得 出 这 样 的 结论 : 在 考虑 随机 品 
声 的 情况 下 , 选择 安培 定律 的 积分 算法 , 并 将 积分 环 路 扩大 到 第 二 个 环 路 积分 时 计算 的 视 向 
电流 密度 分 布 图 更 清晰 ， 而 且 电 流 精细 结构 不 失真 。 由 图 4 可 以 看 出 ， 泡 斑 带 与 电流 带 的 爆 
发 位 置 和 形态 极其 相似 , 所 以 , 计算 的 视 向 电流 密度 分 布 对 于 准确 预测 焰 斑 爆发 的 位 置 和 形 
态 具 有 非常 重要 的 意义 。 为 了 从 理论 上 分 析 积 分 路 径 是 否 合适 ,模拟 时 可 以 在 矢量 磁场 中 设 
定 一 个 电流 密度 的 真实 值 , 但 本 文中 的 数据 是 通过 空间 采样 所 获得 的 观测 数据 , 由 于 受 随 机 
噪声 的 影响 无 法 确定 电流 密度 的 真实 值 ,所 以 文中 是 用 实际 观测 的 数据 验证 ”模拟 方法 的 准 


确 性 ， 即 用 扩大 积分 环 路 (本 文中 为 中 环 路 ) 的 方法 得 到 的 电流 密度 分 布 图 效果 最 好 。 同 时 
我 们 可 以 利用 所 计算 的 视 向 电流 密度 分 布 预 测 泡 广 爆 发 的 可 能 位 置 和 形态 , 也 验证 了 计算 视 


向 


J 电流 密度 的 方法 的 有 效 性 和 可 靠 性 。 未 来 我 们 将 利用 AIA 其 它 波段 的 观测 数据 对 多 个 活动 
区 的 炮 广 带 和 电流 带 进 行 研究 ， 以 检验 该 方法 的 合理 性 和 普 适 性 。 
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Abstract: The explosion of energic solar activity is closely related to the current structure in the 


active region, Ampere's law E (V x B), is the theoretical basis for measuring the current 


density in the active region. Due to the random noise inevitably exist in the measured vector 
magnetic field, so the current density calculated by using the different form of Ampere's law will 
be significantly different. In order to compare the differences between the calculated results of 
different form and explore one of the most practical current calculation methods, this article is 
based on the vector magnetogram in activity region AR11158 which measured by SDO (Solar 
Dynamic Observatory)/HMI (Helioseismic and Magnetic Imager) on February 15, 2011. 
Calculating the distribution of the vertical electric current density in the active region by using 
differential algorithm and integral algorithm of Ampere's law. The results demonstrate that the 
distribution of vertical electric current density j; which obtained by the differential algorithm is 
much more affected by random noise than the result obtained by the integral algorithm, the current 
structures in the current distribution diagram obtained by the integral algorithm is much clearer 
than the former . In addition, the noise signal of the calculated current distribution map will 
decrease sharply, as enlarging the radius of the integral loop, the current structure in the obtained 
current distribution map will become clear. However, when the radius of the integral loop continue 
to expand, part of the fine structures of the current distribution map will also be distorted when 
obtaining a clear current distribution diagram. The result of this study demonstrates that the map 
of current distribution which calculated by expanding the integral loop properly can reduce the 
effect on calculation results from random noise, and obtaining a clear and true vertical electric 
current distribution map. Though the noise is eliminated if the radius of the integral path is too 
large, some fine structures in the current distribution will be lost. Therefore, during the course of 
actual calculating electric current. We should use high-resolution vector magnetogram and use the 
integral algorithm of Ampere's law to calculate the vertical electric current in the active region by 
selecting an appropriate integration path, which can help us to explore the relationship between 


the eruptive flares and the structures of current in the active region. 
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